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ABSTRACT 

Context. Secondary eclipse measurements of transiting extrasolar planets with the Spitzer Space Telescope have yielded several direct 
detections of thermal exoplanet light. Since Spitzer operates at wavelengths longward of 3.6/um, arguably one of the most interesting 
parts of the planet spectrum (from 1 to 3 yum) is inaccessible with this satellite. This region is at the peak of the planet's spectral energy 
distribution and is also the regime where molecular absorption bands can significantly influence the measured emission. 
Aims. So far, 2.2/im K-band secondary eclipse measurements, which are possible from the ground, have not yet lead to secure 
detections. The aim of this paper is to measure the secondary eclipse of the very hot Jupiter TrES-3b in K-band, and in addition to 
observe its transit, to obtain an accurate planet radius in the near infrared. 

Methods. We have used the William Herschell Telescope (WHT) to observe the secondary eclipse, and the United Kingdom Infrared 
Telescope (UKIRT) to observe the transit of TrES-3b. Both observations involved significant defocusing of the telescope, aimed to 
produce high-cadence time series of several thousand frames at high efficiency, with the starlight spread out over many pixels. 
Results. We detect the secondary eclipse of TrES-3b with a depth of -0.241±0.043% (~6cr). This corresponds to a day-side brightness 
temperature of T B (2.2/jm)=2040±185K, which is consistent with current models of the physical properties of this planet's upper 
atmosphere. The centre of the eclipse seems slightly offset from phase <^>=0.5 by A0=-O.OO42±O.OO27, which could indicate that 
the orbit of TrES-3b is non-circular. Analysis of the transit data shows that TrES-3b has a near-infrared radius of 1.338±0.016Rj„ / ,, 
showing no significant deviation from optical measurements. 
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1. Introduction 

During the secondary eclipse of an extrasolar planet (the 
moment the planet moves behind its host star), the contri- 
bution from the direct planet light can be measured. Until 
now such measurements have been the domain of the Spitzer 
Space Telescope, which have led to the first detections of 
thermal emission from ex trasolar planets (Charbo nneau et all 
2005; iDeming et al.l 120051) . Until now a handful of planets 
has been observed at several wavelengths between 3.6//m and 
24yum, with one of the highlights so far being the measure- 
ment of the da y/night-side temper ature contrast of the exoplanet 
HP 189733b (iKnutson et al.ll2007h . 

Unfortunately, a very interesting part of the broad spectral 
energy distribution of hot Jupiters between 1 and 3 //m is short- 
ward of the observing window of Spitzer. This region contains 
the overall peak of the planet's emission spectrum, and can also 
be strongly influenced by molecular bands (e.g. from water, 
methane and CO) in the planets' upper atmosphere. However, 
ground-based observations, which could access this spectral re- 
gion, have thu s far proven very challe nging. Several methods 
have been tried. Richar dson et al.l (l 2003l) used a differential spec- 
troscopic method, searching for the wavelength dependent shape 
of the planet spectrum, but no eclipse was detected. In addition, 
observations of the secondary eclipse of Tr ES- lb in L-band with 
the NIRI spectrograph on Gemini North by lKnutson et al.1 {2007) 
only resulted in an upper limit. 

Photometric attempts to detect the secon dary ec li pse of 
extrasolar planets have been presented by ISnellenl 
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(2005) 



and Sne llen & Covind (120071) . Two partia l eclipses o f 
HD209458b were observed in K-band by ISnellenl d2005h . 
while the telescope was significantly defocused to avoid 
saturation of the array by the bright star. Although milli- 
magnitude precisions were reached, nothing was detected, 
probably beca use of the la ck of baseline on both sides of 
the eclipse (ISnellenl 120051). Subsequ ently, OGLE-TR-113b 
was observed (ISnellen & Covind [2007), which has the great 
advantage that it has several other stars in its field. It was 
shown that randomly offsetting the telescope randomises the 
photometric errors down to 0.1-0.2% per hour. It resulted in a 
tentative detection of its secondary eclipse of 0.17±0.05%. 

In this letter we present our results for K-band photom- 
etry of_fhe_tansn^md_s^condary eclipse of the planet TrES- 
3b dO'Donovan et al.l 120071) . This planet is significantly more 
suitable for secondary eclipse photometry than previous targets, 
since it is in a very close orbit (P~ 1.31 days) and its radius is 
significantly inflated (R~ 1 .3R j up ). Furthermore there is a nearby 
reference star, which can be used for differential photometry. In 
section|2]we present our observations and data reduction. In sec- 
tion|3]we present and discuss our results. 



2. Observations, data reduction and analysis 

2.1. The transit of TrES-3b 

We have observed the transit of the exoplanet TrES-3b with 
the United Kingdom Infrared Telescope (UKIRT) using its Fast 
Track Imager (UFTI: iRoche et alJl2003l) on June 20, 2008. The 
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Fig. 1. K-band lightcurves from the UKIRT transit observations. 
The top panel shows the lightcurve of the reference star and 
TrES-3b before correcting for the positional dependence of the 
flux. The middle panel shows the fluxes of both stars corrected 
for this effect. The dashed line is a polynomial fit to the flux of 
the reference star. The bottom panel shows the final corrected 
unbinned lightcurve of TrES-3b. 

observations, carried out in queue scheduling, lasted for al- 
most 4 hours, starting ~0.6 hour before ingress and ending 
~1.7 hours after egress. The field of view of UFTI is 93x93 
arcsec, with a pixel scale of 0.091 arcsec/pixel. Since the dis- 
tance between TrES-3 (K = 10.61) and the nearby reference star 
(2MASS J175225. 15+373422.1, K-9.77) is ~4 arcminutes, we 
were forced to alternate observations between the target and 
the reference, in order to correct for time dependent atmo- 
spheric and instrumental effects. To reduce overheads from de- 
tector read-out, the detector was windowed to a single quadrant 
(512x512 pixels, 46.6x46.6 arcsec) which increased the cycle 
speed. During each sequence we repeated a nine point dither 
pattern three times before switching to the other star. The ex- 
posure time was 5 seconds for the TrES-3, and 2 seconds for the 
reference star. With a typical overhead of 5 seconds per frame, 
the time required to complete an entire target-reference sequence 
was ~8 minutes. 

The telescope was significantly defocused to 1) reduce in- 
trapixel variations, 2) minimise the influence of flatfield inaccu- 
racies, and 3) prevent the count levels from reaching the non- 
linear range of the detector. Due to the fact that the telescope 
optics had not yet been realigned after the recent switch to the 
Cassegrain focus, the defocused telescope produced a signifi- 
cantly asymmetric PSF, which could be covered by a circular 
aperture of ~40 pixels (~4") in radius. 



Fig. 2. K-band lightcurves from the WHT secondary eclipse ob- 
servations. The top panel shows the raw lightcurve for the refer- 
ence star and TrES-3b. The middle panel shows the lightcurve of 
TrES-3b divided by that of the reference star. The bottom panel 
shows the final, fully corrected, lightcurve. 



The data was flatfielded using a sky fiat constructed from a 
series of frames taken just before the start of the observations of 
TrES-3. No dark subtraction was performed, because any dark- 
current was seen as a contribution to the sky background and 
removed together with the sky. The positions of the hot and cold 
pixels were determined from both the flat field and separate dark 
frames. These pixels were replaced by the values of 3rd order 
polynomial surfaces fitted to the 7x7 pixels surrounding these 
pixels. 

To determine the sky background, first all stars in the field 
were masked using circular masks of 30 pixels. A median sky 
value for the masked image was calculated and subtracted in or- 
der to remove temporal fluctuations in the background. The final 
sky level was determined by averaging the 27 masked frames 
from one cycle, and removed from the images. 

Subsequently, aperture photometry was performed using the 
APER procedure from the IDL Astronomy User's LibraryS with 
an aperture radius of 42 pixels. Any possible residual sky fluctu- 
ations were corrected for by measuring the background value in 
a 42 to 70 pixel annulus around the object. 

The results from the aperture photometry can be found in the 
top panel of Figure Q] Both the flux from the target and the ref- 
erence star was found to be a function of dither position, there- 
fore the data were normalised for each dither position separately. 
These corrected lightcurves are shown in the middle panel of 

1 http://idlastro.gsfc.nasa.gov 



Figure [TJ To remove time dependent effects, we fitted a high or- 
der polynomial to the binned lightcurve of the reference star. 
This fit was used to correct the lightcurve of TrES-3. The final 
unbinned lightcurve for TrES-3 is shown in the bottom panel of 
Figure [TJ 

2.2. The secondary eclipse of TrES-3b 

On July 3-4, 2008 the secondary eclipse of TrES-3 was ob- 
served with the Lo ng-slit Intermediate Reso lution Infrared 
Spectrograph (LIRIS : lAcosta-Pulido et al.ll2002h at the William 
Herschell Telescope (WHT) on La Palma, with the observa- 
tions lasting for the entire night, starting about ~3 hours before, 
and ending ~4 hours after the expected centre of the secondary 
eclipse (at 1:31UT on July 4th). We used the full 1024x1024 
pixel array for these observations, which combined with the 
pixel size of LIRIS of 0.25"/pixel, yields a field of view of 
256x256 arcseconds. This larger field of view enabled us to 
place both TrES-3 and the reference star simultaneously on the 
detector, making these observations ~80% more efficient rela- 
tive to those obtained with UKIRT. We used a 9 point dither pat- 
tern, including small random offsets from the nominal positions 
of each dither point. We used an exposure time of 7.5 seconds, 
with typical overheads of 7 seconds per frame. 

As for the UKIRT observations, we defocused the telescope, 
which resulted in a ringshaped PSF with a radius of about 4 pix- 
els. 

Since the LIRIS detector is non linear at high count levels, 
we applied the following empirical non-linearity correction to 
the data before flat fielding: 
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with F trU e the flux after correction, F meas the measured flux and 
c — 6 ■ 10~ 29 a constant describing the strength of the correction 
term. The value of c was determined by minimising the noise of 
the reference star over the entire night, and comparing this value 
to non-linearity measurements taken by the LIRIS tearrQ. 

All frames were corrected for crosstalk between the quad- 
rants, which was found to be present in a co-add of all the re- 
duced science frames at a level of 10~ 5 of the total flux along 
a row of a quadrant. This correction was done on a row by row 
basis. A flat field was created from a series of domeflats taken 
with the domelight both on and off. These were subtracted from 
each other (to eliminate structure from the emission by the dome 
and telescope) and averaged. 

The small number of hot and cold pixels in the image were 
replaced by a median of the neighbouring pixels in a 3x3 grid. 
The sky was determined from the combination of the 9 images 
in each dither sequence, for which all the stars were masked out. 
We used a circular aperture with a radius of 13 pixels to deter- 
mine the stellar flux and an annulus of 30 to 90 pixels to deter- 
mine the level of possible residual flux in the background after 
sky-subtraction. 

The individual raw lightcurves of TrES-3b and the reference 
star are plotted in the top panel of Figure |2] Before further anal- 
ysis we divided the lightcurve of TrES-3 by that of the reference 
star, shown in the middle panel of Figure[2] The relative fluxes of 
the two stars are correlated with their positions on the detector. 
These effects are corrected for by fitting linear functions between 
the x,y-positions and the relative flux for data outside of the ex- 
pected eclipse. These fits are then applied to the whole dataset. 




Fig. 3. The 9-point binned lightcurve of the transit of 
TrES-3b. The solid line shows the best fitting model with 
R/>=1.338±0.016R/ U/ ,. The bottom panel shows the residuals af- 
ter subtracting this model fit. 
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Fig. 4. Final, 9-point binned lightcurve of the secondary eclipse. 
The solid line is the best fitting model with AF=-0.241±0.043% 
and an offset from phase 0.5 of A0=-O.OO41+O.OO18. The bot- 
tom panel shows the residuals after subtracting the best fitting 
model. 



In a similar way a small dependence on airmass was removed. 
The corrected unbinned lightcurve of TrES-3b can be found in 
the lower panel of Figure [2] 



3. Results and discussion 

3.1. The transit of TrES-3b 

The final transit lightcurve binned per 9-point dither cycle is 
shown in Figure[3] The RMS-noise level out of transit is 0.0060, 
while the theoretical noise level for these observations is 0.0029, 



which is actually dominated by the detector's readnoise of 
26e~/pixel, due to the large aperture used. 

Since a dependence of the transit depth on wavelength can 
reveal characteristics of the planet atmosphere, we fitted the ra- 
dius of TrES-3b, keeping the parameters of the host star and 
the planet's orbit (mean stellar density, orbital i nclination and 
ephem eris) fixed with respect to the values of ISozzetti et al. 
(2008). For the limb darkening we used parameters from lClareti 
d2000 ). assuming a temperature of T — 5750K for the host star of 
TrES-3b, a metalicity of log[Fe/H]= -0.2 and a surface gravity 
of log(g)= 4. 5. The transit w a s mod elled using the IDL pro- 
cedures from Man del & Ago! j2002), resulting in a radius of 
1.338±0.013Ry„ p at ^ 2 /v =1.58. The reduced^ 2 is significantly 
larger than unity. Since we fitted to the binned lightcurve, this in- 
dicates that the noise does not scale with yfN, meaning that we 
have an additional, unknown, source of correlated noise. To ob- 
tain a better estimate of the uncertainty in the planet radius, we 
forced x 2 / v = 1 by scaling up our errors in the binned lightcurve 
uniformly. In this way we obtain Rp=1.338±0.016Ry„ p . The K- 
band radius shows no significant deviation from the op t ical ra - 
dius of 1.336^iR /ap , as measured by ISozzetti et all (120081) . 
The uncertainty in our K-band radius appears smaller, but this is 
solely due to the fact that we keep the parameters for the host 
star and the planet's orbit fixed. 

We show that we can obtain the K-band radius of an exo- 
planet at a precision of This could be further improved by 
observing the reference star and target at the same time both be- 
cause of the increase in the number of frames during transit and 
because rapid atmospheric fluctuations can be corrected for. 

3.2. The secondary eclipse of TrES-3b 

The fully corrected lightcurve binned by one 9 point dither cy- 
cle is shown in Figure [4] Before binning, and after subtracting 
the best fitting model (see below), we clipped the lightcurve at 
±0.01, removing 43 of the 1800 points (2.3%), of which 3 points 
(0.9 %) during the ec l ipse. T he binned lightcurve was fitted using 
the iMandel & Agoll (120021) procedures, with all the parameters 
fixed to the transit model, except the centre and the depth of the 
eclipse (of course without limb darkening). 

The best fitted lightcurve is over-plotted in the top panel of 
Figure 2J with the bottom panel showing the residuals. The out 
of eclipse RMS noise is 0.0038, while the expectation from the- 
oretical noise statistics is 0.0015. The dominant source of noise 
in these noise calculations is the sky contribution. The total flux 
from the sky within the aperture is, even for the brighter refer- 
ence star, ~50% higher than the stellar flux. 

The fitted depth of the secondary eclipse is 
AF =-0.241+0.021%, with a reduced chi-squared of 
X 2 /v=1.72. As for the transit data, the reduced x 2 f° r the 
fit to the binned lightcurve is greater than unity, indicating that 
we have some extra systematic noise. We also find a timing 
difference in phase of A0=-O.OO41±O.OO13 (at ^ 2 /v=1.72). The 
offset in the timing is intriguing, because it may indicate that 
the orbit of TrES-3b is not circular but slightly eccentric, with 
ecosw =-0.0066±0.0021. Although the formal significance of 
the offset is ~3cr, the residual systematic noise in the lightcurve 
decreases its significance, in particular due to several points 
around <p = 0.49. We also refitted the data with the centre of 
the eclipse fixed to 0=0.5, and obtained in this case a depth of 
-0.233±0.020 at* 2 /v=1.77. 

To test the influence of the more deviated points around 
0=0.49, we repeated the fit while excluding them. The depth 



subsequently becomes -0.190^[H51 with an offset in phase of 
A(p =-0.0036+™° 22 Otf 2 /v=l. 39). Although the significance of 
the fit is reduced, we still detect the secondary eclipse at a 6<x- 
level. 

Because it is clear that red noise dominates the uncertain- 
ties in our measurements, we tried to characterise the uncer- 
taint y using the "resi dual permutation" bootstrap method as used 
by IWinn et al.1 |2008). The best fitting model was subtracted 
from the data, after which the residuals were shifted between 
1 and 1800 points and added to the model light curve. These 
light curves were refitted for each individual shift. The result- 
ing distributions of parameter values represent better estimates 
of the real uncertainties in the data. With this analysis we find an 
uncertainty of 4.3- 10~ 4 in the eclipse depth and 2.7- 10~ 3 for the 
offset in phase, which we use in the remainder of the paper. 

Assuming a blackbody model for the thermal emission spec- 
trum of TrES-3b, we can estimate the day-side brightness tem- 
perature in the K-band. For the ste llar temperature, we use 
T eff =5650±75 K dSozzetti et al.l2008l) . From the fit to the eclipse 
depth using the full lightcurve we obtain T B =2040±185K (90% 
confidence interval). 

Due to the high incident stellar flux, For tnev et alj (2008) 
classify this planet as belonging to the pM class of planets. This 
class is expected to have an inversion layer in their upper at- 
mospheres and to show a large day/night contrast. Comparing 
the depth of the secondary eclipse with their model (their Figure 
14), we see that our results are indeed consistent with this model. 
However, since we have measured the brightness temperature 
of TrES-3b at only single wavelength, Spitzer measurements at 
longer wavelengths are needed to confirm the presence of an in- 
version layer and thus classify TrES-3b as a pM planet. 

4. Conclusions 

We have detected the secondary eclipse of the exoplanet TrES- 
3b in K-band at a level of -0.241 ±0.043%, indicating a bright- 
ness temperature of TB(2.2yum)=2040±185K. The eclipse tim- 
ing shows a small offset from 0=0.5, which would imply a non- 
circular orbit, but this needs to be confirmed by Spitzer observa- 
tions. Also the radius of TrES-3b is determined in K-band to be 
1.338±0.016Ry l(p , showing no deviations from optical measure- 
ments. 
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